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Abstract 

We study the matter effects to neutrino oscillations — * v e and Vn — * v„ at very long baselines. 
It has been observed that, for a very long baseline L ~ 10000 km, the resonance peak for the 
Vu, v e oscillation and the local minimum of the v„ survival probability occurs at similar energies. 
With a good knowledge on the absolute value of Am^, measurements of the above oscillation 
probabilities, P^ e and Pau, can be performed with the neutrino energy tuned to the resonance 
peak of the Va — > v e oscillations. We show that the variations of CP violating phase and the solar 
neutrino mixing parameters have negligible effects on P^ e and Pun for such a long baseline. Hence 
P^ e and Pan together determine the mixing angles #13 and #23- Around the resonance peak for 
v a v e oscillations, the dependencies of P^ e and P^ on the above mixing angles are worked out 
in detail, and the implications of such dependencies are discussed. 
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The understanding of neutrino masses and mixing matrix is crucial to unveil the mystery 
of lepton flavor structure. The updated SKK analysis of the atmospheric neutrino data gives 



1.9 ■ 10~ 3 eV 2 < | Am^ | < 3.0 • 1(T 3 eV 2 , sin 2 29 23 > 0.9. (1) 

This is a 90% CL. range with the best fit values given by sin 2 2^3 = 1 and Am| x = 
2.4 • 10 -3 eV 2 respectively. The scenario of — > v T oscillation for atmospheric neutri- 
nos has been confirmed by the K2K experiment |2|. Furthermore the results in the solar 
neutrino oscillation measurements are also confirmed by KamLAND reactor measurements 
. Combining these measurements, the LMA solution of the solar neutrino problem is 
established and the updated 3a parameter ranges are given by P 

7.4 • 10~ 5 eV 2 < Ar/4 < 9.2 ■ 10~ 5 eV 2 , 0.28 < tan 2 9 12 < 0.58, (2) 

with the best fit values Am 2 x = 8.2 ■ 10" 5 eV 2 and tan 2 12 = 0.39. 

Despite the achievements so far in measuring the neutrino mixing parameters, the sign 
of Am| 1; the mixing angle #13 and the CP violating parameter Sqp in the mixing matrix 
remain to be determined. The sign of Am^ can be determined through matter effects to 
Vn — > v g {pp — > v e ) and v e — > {p e — > z/ M ) oscillations in the very long baseline experiments 
la, Ml. However, such oscillations are sensitive to the mixing angle #13, which has been 
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constrained by the reactor experiments (8J, |9J. The CHOOZ experiment [8J] gives a more 
stringent constraint on # 13 with sin 2 26* 13 < 0.1 for a large Am^ (90% C.L.). A recent global 
fit gives sin 2 2^i3 < 0.09 (0.18) at 90% CL. (3a) Q- 

It is known that the oscillation — > z/ e , for example, is enhanced (suppressed) by the 
matter effect for Am 2 x > (Am^ < 0). On the contrary, the — > v e oscillation is enhanced 
(suppressed) by the matter effect for Am^ < (Am^ > 0). In the following discussions 
we shall assume Am^ > unless explicitly specified. It is noteworthy that, although the 
matter enhanced angle 8^ reaches to 7r/4 at the resonance energy, the oscillation probability, 
Pivp, — > v e ) = Pn&j remains to be small unless for a very long baseline. Indeed, it is 
demonstrated that [llj, in the constant density approximation for the Earth density profile, 
the condition for a maximal P Me reads: 

pL max _ (2W + 1)7T y g lg 1Q 3 km . g / cm 3 ; (3) 

tan 20i3 

where p is the averaged mass density of the Earth, L max is the required oscillation length, 
and n is an integer. For n = and sin 2 2# 13 = 0.1, L max is found to be 10200 km [12]. A 



smaller sin 2 29\z requires an even larger oscillation length. Clearly, to take a full advantage 
of the matter enhancement, one is led to consider experiments with a baseline of the order 
of the Earth diameter 

It has been demonstrated that the matter effects to P{Pn — > v T ) = P,,.t and P(vu — > 



Ufj) = P^/j, can also be significant in the vicinity of P^s resonance peak 
the criterion for the maximal matter effect to P Mr is identified to be E res : 



QQ 



. In fact, 



' -^pcak where the 

former is the resonance energy for the oscillation P„ e while the latter is the energy that the 
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vacuum approximation of P^ T attains its peak value. Such a criterion amounts to 

pL max ~ (2n + l)7r(cos 20i 3 ) x 5.18 • 10 3 km • g/cm 3 . (4) 
0.1, L max ~ 4400 km, while, for n = 1 with the same sin 2 29 
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9700 km. In the latter decreased from its vacuum oscillation 



For n = and sin 26* 13 
value, L max 

value by as much as 70%. This implies that the muon neutrino survival probability P^ also 
receives significant matter effects in the same energy range. 

It is interesting to note that the matter effect to P^ and that to P^ are both significant in 
the same energy range for baselines of the order 10 4 km. Hence simultaneous measurements 
of P^ and P^ in such a baseline provide a stringent test on the matter effect, which is the 
focus of this work. To discuss these oscillations, we begin with the relation connecting flavor 
and mass eigenstates of neutrinos, v a = YliUai^i, with U the Pontecorvo-Maki-Nakagawa- 
Sakata mixing matrix given by 



u 



C12C13 



S12C13 



-S12C23 - Ci 2 Si 3 S23e i,5cP C12C23 - Si2Si 3 S2 3 e i ' 5cp 



(5) 



si 3 e 
C13S23 

S12S23 — C12S 130236'"^ — C12S23 — Sl2 s 13 c 23 e ''"''~ J± ' C13C23 

where s i3 - and Cy denote sin^- and cos 6^, respectively. For the Dirac type CP-phase Sep 
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o*<5cp 
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we allow its value ranging from to 2tt. The evolutions of neutrino flavor eigenstates are 
governed by the equation 
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where \v(t)} = (u e (t), u^t), u T (t)) T , Am. 



m| — m 2 is the mass-squared difference between 
the i-th and j-th mass eigenstates, and V = \^2GfN £ is the effective potential arising 
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FIG. 1: The Am^ dependence of the transition probability for sin 2 2# 13 = 0.1, sin 2 26*23 = 1) 
5cp = 0, and the baseline length L = 9300 km. 

from the charged current interaction between u e and electrons in the medium with N e the 
electron number density. Numerically V = 7.56 x 10~ 14 (p/[g/cm 3 ])Y" e [eV] with Y e denoting 
the number of electrons per nucleon. For the Earth matter, Y e ~ 0.5. One solves Eq. (JBJ) 
by diagonalizing the Hamiltonian on its right hand side. This amounts to writing the right 
hand side of Eq. © as U' H'U'^\v(t)) with U' the neutrino mixing matrix in the matter and 
H' = diag(J5i, E 2 , E 3 ) the Hamiltonian after diagonalization. To obtain various oscillation 
probabilities described later, we have used the parametrization in [15J for the Earth density 
profile [16j . For illustration, we shall take L = 9300 km, which is a distance between Fermilab 

n 

and Kamioka [13J]. This distance is in fact very ideal for studying P Me and P MAt since both 
oscillations receive significant matter effects around such a distance. 

To measure and near the resonance peak of the — > u e oscillation, it is crucial 
to know the precise value of Am^. The current range for this parameter, Eq. ((TJ, results 
in a corresponding range for the peak energy of P Me as shown in Fig. ^ The energy value 
at the peak of P^ e is proportional to Am^ as dictated by the resonance condition 2E U V = 
Am 2 ! cos #i3 while the maximal magnitude of P^ e is independent of Am^. This is e asih 
understood by inspecting the expression for P^ e in the constant density approximation 



asily 



4 



E [GeV] 



FIG. 2: The transition probabilities P^ e and P w for a baseline length L = 9300 km with sin 2 29\ 3 = 
0.1, sin 2 2^23 = 1 and 5 cp = 0. The gray zone represents the region with the neutrino energy E v 
ranging from 5.5 GeV to 6.5 GeV, where the maximum of P^ e and the local minimum of Pan are 
located. 



P Me = sin 2 # 23 sin 2 20™ sin 2 (1.27 A™ L/E u ) , (7) 

with L in km, E u in eV, 



A™ = y (Am§! cos 29 13 - 2E u Vy + (Am 2 31 sin 29 13 ), (8) 

and 

_ sin 20i3- Am 2 , 

sm20 13 - — . (9 J 

^31 

Clearly, at the resonance energy, the last factor on the right hand side of Eq. (JJJ) does not 
depend on Am^. 

As pointed out earlier, the matter effects to P^ e and P w are both significant in the same 
energy range for our interested baseline length. This is indicated by the gray zone of Fig. El 
where the energy dependencies of P^ e and ? w are depicted. The probability P Ufl shows an 
intriguing feature in this energy zone. Due to the large matter effect, P MAt does not drop to 
zero while reaching to its local minimum. For Am 2 , = 2.4 x 10~ 3 eV 2 , sin 2 2# 13 = 0.1 and 




FIG. 3: The contour graph of the transition probabilities P ae (solid line) and Pan (dashed line) in 
the sin 2#i3-cos 2023 plane for <5cp = and the baseline length L = 9300 km. The plotted P^ and 
Pau are the averaged probabilities for neutrino energies ranging from 5.5 GeV to 6.5 GeV 

sin 2 2#23 = 1, the probability P ue peaks at 6.4 GeV while the local minimum for P au is at 
6.1 GeV. We have varied the values for #13 and #23- The energy value for the maximal P ue 
and that for the local minimum of P uu are shifted. However, the shifts are within 300 MeV 
for both cases. 

The probabilities P ua and P ae as functions of #23 and #13 are shown in Fig. El The range 
for sin 2^13 is taken from to 0.32 in accordance with the reactor bound sin 2 2^13 < 0.1 
while that for cos2# 2 3 is from —0.32 to 0.32 which is equivalent to sin 2 26* 2 3 > 0.9 [1]. The 
oscillation probabilities plotted here are the averaged ones for neutrino energies ranging from 
5.5 GeV to 6.5 GeV. For Am^ = 2.4 x 10~ 3 eV 2 , both P ue and P au receive large matter 
effects. We argue that these probabilities do reflect the relative event rates. In fact, to 
obtain neutrino event rates, one has to convolve the oscillation probabilities with the initial 
muon neutrino flux and the conversion rates from neutrinos to charged leptons. However, 
since our interested energy range is narrow enough, the neutrino-nucleon scattering cross 
section can be treated as a constant. The same is true for the u u flux. We note that P ue is 
essentially only sensitive to sin26 , 13 for P^ e < 0.1. However P ue is also sensitive to cos2# 2 3 
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FIG. 4: The CP-phase dependencies of the transition probabilities P^ e and P w for sin 2 2#i 3 = 0.1, 
sin 2 2#23 = 1) and the baseline length L = 9300 km. 

as it gets larger. The contours for constant P^'s show more complicated structures. We 
observe that, for a fixed 6 I 13 , a given P m yields two solutions for #23 for P m < 0.15. The 
measurement of P^ e might help to remove this degeneracy if P^ e is large and consequently 
sensitive to cos 26*23. We also note that the above degeneracy no longer exists for P^ > 0.15. 

We point out that simultaneous measurements of P^ e and P MAt determine 6*13 and #23- It 
is found that the variations of Am^ and 812 give negligible effects on the above oscillation 
probabilities. The effects of CP violating phase 5qp on P Me and P MAl are also studied. With 
AfBjj = 8.2 ■ 10" 5 eV 2 , sin 2 26*23 = 1> sin 2 2^ 13 = 0.1, the dependencies of the peak value of 
P Me and the corresponding minimal value of P w on #cp are shown in Fig. 0] It is seen that 
the peak value of P Me differs by no more than 4% for the entire range of #cp. Similarly, the 
corresponding minimal value of P w differs by no more than 3%. 

The matter effect to P^ can be identified by the resonance enhancement to the oscillation 
v \i ~~ * v ei given our assumption that Am^ > 0. On the other hand, the matter effect to P w 
is seen by comparing Fig. El with a corresponding plot, Fig. where matter effects to P MAt 
are neglected. In Fig. EJ one can see that P^ is not sensitive to sin 26 \% and the contours for 
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FIG. 5: The contour graph of the transition probabilities P^ e (solid line) and P^(dashed line) for 
^CP = and the baseline length L = 9300 km. Vacuum oscillation is assumed for P/m. 

constant P^s behave very differently from those in Fig. 01 particularly for a larger sin 28 is. 
We observe that, in the vacuum oscillation case, P MAt < 0.15 for almost all the parameter 
space in the sin 26*i3-cos 2#23 plane. Hence an observation of P M/J significantly greater than 
0.15 is a signature of the matter effect. The parameter space for this to occur is situated in 
the upper right corner of sin 2# 13 -cos 2# 2 3 plane. This is a region with 823 < 7r/4. 

The matter effect to P MM can also be tested in other parameter regions. This is true 
when values of #13 and 823 are accurately measured in other future experiments. There 



are extensive discussions on measuring 813 in reactor neutrino experiments |1 



19j. The 



issue o 
in 
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determining 823 in the three-flavor oscillation framework has also been discussed 
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22j. It is expected that sin 2 26*23 can be measured to an 1% accuracy in the 
J-PARC— ► Super-Kamiokande (JPARC-SK) experiment j^J. This translates into a 10% 
accuracy for the determination of sin 2 823 for 823 around 7r/4 at 90% C.L., provided sin 2 26*13 
is accurately measured by the reactor experiment j^j], say around 2% at 90% C.L. ac- 
cording to the analysis in Q. Given the future accuracies for determining sin 2 6^3 and 
sin 2 26^3, it is possible to identify the matter effect in the muon survival probability P^. 
For illustration, let us assume, for example, P Me and P w are measured to be 0.4 and 0.1 re- 
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spectively in the very long baseline experiment, with some uncertainties understood for each 
measurement. Taking into account the matter effect to P^, these oscillation probabilities 
correspond to (sin 2# 13 , cos 2^23) = (0.26,-0.03) according to Fig. El However, by assum- 
ing vacuum oscillations for P MAt , one obtains (sin 2# 13 , cos 2# 23 ) = (0.22, —0.22) according to 
Fig. |5J Therefore the extracted central values of sin 26 13 differ by 16% with and without the 
matter effects to P MAt taken into account. This discrepancy is considerably larger than the 
expected uncertainty in the determination of sin 2#i 3 by the reactor experiment, which is 
approximately 4% with the above-mentioned uncertainty 5 (sin 2 2# 13 ) ~ 2% and the central 
value sin2#i 3 = 0.26. Hence the matter effect to P MM can in principle be identified by the 
extracted value of sin2# 13 , provided there are sufficient accuracies in the measurements of 
P^e and Pnn- Furthermore, we note that sin 2 # 23 = 0.52 with the matter effect to P M/J taken 
into account while sin 2 623 = 0.61 without the matter effect. The difference on the extracted 
value of sin 2 623 is about 17%, which is also larger than the expected 10% uncertainty in de- 
termining sin 2 623 by the JPARC-SK experiment. Therefore, if one can accurately measure 
P Me and P M(J , the matter effects to P w can also be identified through the extracted value of 
sin 2 023- 

In summary, we have discussed the matter effects to —>■ u e , oscillations at very long 
baselines. We illustrate these effects with a baseline L = 9300 km, the distance between 
Fermilab and Kamioka. Around this baseline length, the matter effects to P Me and P MM 
are both significant near the resonance energy for z/ M — > u e oscillations. We have worked 
out in detail the dependencies of these probabilities on the mixing angles 813 and #23 with 
Am 31 taken to be 2.4 x 10~ 3 eV 2 . We have also shown that these probabilities are not 
sensitive to the variations of other neutrino oscillation parameters. We have identified the 
parameter space in the sin 2#i 3 -cos 2#23 plane in which the matter effect to P MM can be most 
easily identified. We also argued that it is possible to identify the same matter effect in 
the remaining parameter space provided 6*i 3 and 6*23 are accurately measured in other future 
neutrino experiments. 
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